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Abstract 

In this paper we find the inverse and direct recursion operator for 
the intrinsic generalized sine-Gordon equation in any number n > 2 
of independent variables. Among the flows generated by the direct 
operator we identify a higher-dimensional analogue of the pmKdV 
equation. 
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Introduction. General recursion operators as introduced by 
Olver [23] are pseudodifferential operators J2t=- r fi-^x ^* mapping 
symmetries to symmetries, thus capable of generating infinite series of 
them. An important generalization by Guthrie ^3] eliminates inherent 
problems of interpretation of D~ . Moreover, a link exists between 
Guthrie's definition and the theory of coverings (Krasil'shchik and 
Vinogradov |17| I18j ) which implies that main constituents of Guthrie 
recursion operators trivialize in higher dimensions unless the system 
in question is (maximally) overdetermined |20| . 

The goal of this paper is to demonstrate that there do exist overde- 
termined systems in any dimension which admit a nontrivial recur- 
sion operator under Guthrie's definition, thus escaping the bi-local 
approach long established in multidimension ^2 ^] . 

The system in question is the intrinsic generalized sine-Gordon 
equation (IGSGE) [Ij, describing metrics of re-dimensional Rieman- 
nian spaces M of constant sectional curvature K immersed into a 
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(2n — 1) -dimensional Riemannian space M of constant sectional cur- 
vature K. The subcase K = of flat metrics on M is usually re- 
ferred to as the intrinsic generalized wave equation (IGWE). These 
systems are known to be integrable in the sense of soliton theory. 
In particular, these equations posses a Backlund transformation [I] 
as well as soliton solutions [H3 EH- Closely related are the gener- 
alized sine-Gordon equation (GSGE) and generalized wave equation 
(GWE) [Tl 1291 05] . which describe the immersions themselves, as well 
as the Generalized Equation (GE) and the Intrinsic Generalized Equa- 
tion (IGE) jHlEl; referring to pseudo- Riemannian metrics. All these 
systems are integrable as well, with ever growing bibliography of re- 
sults, e.g., dUnillllEU, to list just few. 

Lie symmetries of IGSGE and IGWE were computed by Tenenblat 
and Winternitz [27] along with the corresponding invariant solutions, 
this result being extended to GE by Ferreira However, apparently 
no higher symmetries have been written yet. We utilize the recursion 
operator found in this paper to generate a local flow of third order, 
which we call the generalized pmKdV equation, and conjecture it to 
be an integrable evolution system. 

Unless otherwise stated, indices range between 1 and n. 

The equation. (j3J [201 I27J) Let an n-dimensional Riemannian 
manifold M of constant sectional curvature K be immersed in an 
(2n — l)-dimensional Riemannian manifold M of constant sectional 
curvature K. Then M admits orthogonal coordinates x 1 ,...,^", 
I = (v 1 dx 1 ) 2 + • • • + (v n dx n ) 2 , with the n Lame coefficients v , . . . , v n 
satisfying an algebraic constraint 

i 

and n 3 — n 2 + n first-order differential equations 
v) = v j h j \ j ■■ i. 

h*f = -hf - KvV - h si h sj , i > j, 

h v = - h f - £ ,rh>\ z < j, 

Here i 7^ j, are n 2 — n auxiliary functions (rotation coefficients). 
The lower indices denote differentiation, hence v l j stands for dv l /dx 3 
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etc. 

This system is referred to as IGSGE when K ^ and IGWE when 
K = 0. Our results, formulated for IGSGE, are easily transferable to 
IGWE. In the sequel we apply the convention h w = to remove the 
restrictions on summation indices. 

To start with, we check the highly overdetermined system © for 
consistency. The system is orthonomic in the sense of the Riquier 
theory (see (HI § 6] for a compact summary) under a suitable ranking 
of derivatives. Namely, if derivatives cT lH hrn / '{dx l ) ri . . . (dx n ) rn are 
ranked by their total order r\ + • • • + r n , with ties broken by the degree 
r n in x n , with ties broken by the degree r n _i in x n_1 , etc., then (J2J) 
is (irrespectively of the ranking of the dependent variables v 1 ,^) 
resolved with respect to the highest rank derivatives and hence so is 
its infinite prolongation resulting from differentiation of each equation 
of (j2j) with respect to an arbitrary combination of x*'s. 

Derivatives on the left-hand side of (J2J) and its infinite prolongation 
are called principal and become functions of the remaining - paramet- 
ric - derivatives. Now, the system (j2J) is easily seen to be passive, 
meaning that all its integrability conditions (equality of cross deriva- 
tives) are identically satisfied and are so even without reference to the 
algebraic constraint JQ). Moreover, one easily checks that the deriva- 
tive Yli v%v \ = of Eq. (Q) with respect to an arbitrary x k is satisfied 
by virtue of the first two equations of system @. We then conclude 
that the whole system ((2j) is passive. 

The constraint Q means that (y , . . . ,v n ) lies in 5 n_1 , the unit 
sphere. Then v n can be expressed through v , . . . , v n ~ l on each of the 
two hemispheres v n > and v n < of S n . Hence the complete list 
of parametric derivatives to be used below: 

h lj , i + h 

U ^J 1 > j ' ( 3 ) 

r 
r 

the number r of lower i's or j's being arbitrary positive. 

In terms of geometric theory of systems of PDE 16., we have the 
following picture of the system (J2J). Let Y denote the product 
of the unit sphere 5* n_1 , coordinatized by v 1 , . . . , v n , and the arith- 
metic space coordinatized by , i ^ j. Consider the trivial 
bundle it : Y x M — > M and its infinite prolongation J°°7r. Then 
system (j2J) and its infinite prolongation (consisting of all differential 
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consequences) above determine a submanifold £ in J°°7r. The system 
being passive, differential consequences with identical principal deriva- 
tives on the left-hand side are equivalent. The parametric derivatives 
© then serve as coordinates along the fibres of £ — > M. 

Given a function / G C^J^ir, its restriction f\g to £ is com- 
puted if every principal derivative / depends on is substituted by the 
right-hand side of the corresponding equation (j2J or its differential 
consequence until all principal derivatives disappear (which is guar- 
anteed by the use of Riquier ranking). This leaves us with a function 
independent of derivatives of v l . Then v n are replaced with their 
expression through the coordinates v , . . . ) v n . This is a workable 
procedure that leaves us with a function involving parametric deriva- 
tives (JHJ only. 

The usual total derivatives with respect to x k admit restriction on 
£, namely 

i=l i<j 

■ 1 *'i 'jjjT, +EE /i }..j >■ '77777 ' 

r>0 i>j r >0 i<j 3—3 

r r 

where coefficients h?\s, etc. are determined from system @ and 
its prolongations. 

The zero-curvature representation of system (|2j), due to 

Beals and Tenenblat 4 , consists of sparse antisymmetric 2n x 2n 
matrices A/m, A; = 1, . . . , n, given by 

^ = = {t z ~ z ) ^ ~ Tz 6 ^ (4) 

A n+i,n+j = ^ij _ ^ji 

Here z is the so-called spectral parameter. The zero curvature condi- 
tion A/ k \ x i — An\ x k + [A(fc), A (j)] = 0, which holds as a consequence 
of system (J2J) , is geometrically expressed as 

(A A (fc) - D fc A (J) + [A (fe) , A (l) ])| £ = 0. (5) 

Symmetries and recursion operators. As is well known (see, 
e.g., jlfij). a symmetry of PDE can be identified with a vertical vector 
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field Q = ^2 L Q L d / 'dq L on the corresponding manifold £ — > M. The 
field Q is required to commute with the total derivatives Dk restricted 
to £ . It turns out that symmetries of a system of equations {F l = 0} 
have to satisfy the linearized system {£ F i = 0}, where 

with q l running through the coordinates on the manifold £ (paramet- 
ric derivatives). Moreover, the condition of commutativity with 
implies that the coefficient Q K corresponding to a derivative q K = 
of another coordinate q L is expressible as Q K = D^Q 1 ■ Thus, the only 
coefficients to be found are the Q 3 corresponding to variables q J of 
zeroth order (v % and h 13 in our case). 

Hence, regarding our system (pQ), (J2J), symmetries are determined 
by functions V 1 , H %3 on the manifold £ (i.e., depending on (J5J)) that 
satisfy 

i 

Vf = v 3 H ji + h ji V j , j^i, 
V? = ~Y,(v s H is + h is V s ), 

//'' = -Hf - KirV + viV 4 ) (6) 
- (h si H s i + h S3 H si ), i > j, 

Hf = -Hf - ^ (h is H js + h js H is ), i < j, 

Hl j = h ik H kj + h kj H ik , j^i^k^j, 

where VJ = DjV 1 etc. 

Following Guthrie ^3] and our earlier observation 20, , we interpret 
recursion operators as Backlund autotransformations for the linearized 
system (JJJJ). 

As a rule |23[ [2j, there exists a recursion operator that can be 
written in terms of an auxiliary system of equations closely related to 
the zero-curvature representation, namely 

W^ = [A W ,W]+4 W W, (7) 

with ^A( fe ) computed componentwise. In our case, W is an antisym- 
metric 2n x 2n matrix, akin to matrices Aa.). Compatibility of system 
(|7jl follows from the zero curvature condition 
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Now, it is easily checked that if V l ^H %3 are symmetries and W 
satisfies (0, then V'\H' ij given by 



V H = 2zV i - 2z^2v s W is , 

(8) 

H*j = Kv j v s W s ' n+i - {\K + 2z 2 )W j ' n+i 



are symmetries as well. Thus, formulas (0 and (JBJ) determine a family 
of recursion operators 1Z Z for the IGSG equation. When applied to a 
(local) symmetry of IGSG, 1Z Z produces as much as \n{n — 1) symme- 
tries, all of them highly nonlocal as a rule. In such cases the conven- 
tional recursion operator is to be found among inverses of 1Z Z — Aid 
for suitable A and appropriate choice of the spectral parameter z. 

The conventional recursion operator turns out to coincide 
with TZ~ l whenever K ^ ±Az 2 . Inversion of 1Z Z (for arbitrary n) 
uses , (jHJ , the linearized equation (jHJ) and the same linearized equa- 
tion © written for V h ,H' %:> (referred to as (JHI)) to express V l ,H 13 in 
terms of V H ,H H:) while eliminating as many W %3 as possible. 
To start with, we find 



2z ^ 



from the first equation of ((HJ) and 

H ij = -DiW ij + h is W js 

s 

+ (r~ ~ t"^) 2 " z ) wj,n+i + ^-rf^w^ 

from the upper left off-diagonal part of 0. For all i ^ j, we then 
have 

W i,n +j = 2 H'» + AKv* J2 S v S H" JS 



K + Az 2 (K + Az 2 )(K -Az 2 -2K{vi) 2 ) 
2Kv i v j W j ' n+j 



K - Az 2 -2K(vi) 2 



from the second equation of (|HJ), assuming that K ^ —Az 2 . Simi- 
larly, all remaining components W %3 and W n+l ' n+:! can be expressed 
algebraically in terms of W l ' n+l , assuming additionally that K ^ Az 2 . 
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Omitting details of tedious computation, we finally obtain W l ' n+l in 
terms of "potentials" Q l , which can be characterized as follows. Let 

q i = K(v i ) 2 + Y((h si ) 2 + (h* s ) 2 ), 

s 

q\ = -2h ik h ki for i ^ k. 

Then Diq l k = D^q\ for all k ^ I on £, meaning that horizontal forms 
£ J = ql dx k determine an abelian covering over (see ^S] ) . Lin- 
earized forms = Ql dx k with 

Qj = 2Kv i V' i + 2 s P{h si H' si + h is H' is ), 

(9) 

Ql = -2h lk H' kl - 2h kl H nk for i ^ & 

then satisfy DiQ l k = D^Q\ on (JJ3), hence determine an abelian cover- 
ing over the linearized equation (0). The required potentials Q % are 
just naturally chosen nonlocal variables of this covering and may be 
thought of as independent quantities satisfying 

D k Q l = Ql (10) 
In terms of potentials Q l , the final result of inversion of 7Z Z is 
V i = 2^ v s D % H' is + 2^h ti Y^ v s H' ts ~2J2 hti Yl 

s is is 

- 2K(v i fV H + ^2v s (h is Q i - h si Q s ), 

s 

H ij = -2D 2 H Hj -2^2h is D s H' sj - 2^K{v i ) 2 + ^{h^f^H^ 

s s 

-2^2^h is H' is + 2h ij h si H' si + (hf + h u h st ^H' s i 
-2Kh ij (2v i V' i -v j V' j ) 

- h^Q 1 + (hf + hiShjs ) Q j - hiShSj Q s 

s s 

(11) 

(terms containing z divided out). Equations (0), (|9|)- (fTT|) imply @. 
That is, (f9*])-()ll|) determine a recursion operator, henceforth denoted 
C, for the IGSG equation. 

Remark 1 The total derivatives in Q1UJI are those of the covering 
equation (they commute by virtue of DiQ l k = D^Q]) and hence the 
equation (jTU|) is somewhat formal. Only when C is to be applied to a 
seed symmetry a with components V ll ,H l13 , then the Q k , as well as 
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= Q l k dx k , induce objects on £ (or an appropriate covering thereof 
if the seed is a nonlocal symmetry). It happens rather often that 
the induced forms are exact on £, hence the potentials Q l exist as 
functions on £ and the resulting symmetry C(a) is local. If are 
not exact on £, then they induce an abelian covering £ over £ and 
C{<r) is then a nonlocal symmetry with shadow defined on £ |18l I2(J| . 
(This usually happens with 1Z z (a) and explains why operators 1Z Z fail 
to generate local symmetries.) 

Other cases. Routine computation shows that 7Z Z — Aid is invert- 
ible for all A 7^ 0, but symmetry-generating properties of (1Z Z — A Id) -1 
are no better than those of 7Z Z itself. Therefore we continue with 
A = 0. 

In two singular cases, when K = ±Az 2 , the Backlund autotrans- 
formation 1Z Z admits a reduction of the set of nonlocal variables (see 
|21j for generalities). If K = — 4z 2 , then 1Z Z can be written in terms 
of 2re nonlocal variables p % = Y^j=i v^W 13 , i = 1, . . . ,2n as 

V* = 2z(V i - p*), H' ij = Az 2 v j p n+ \ 

with p % being subject to another determining linear system p l xk = 
B^p 3 + B'^y whose exact form is not important here. It follows 
immediately from the second equation that H Hl /v 1 = • • • = H' nl /v n ; 
therefore the image of TZ Z is not the full £ , hence 1Z Z is noninvertible. 

When K = Az 2 , a maximally reduced set of nonlocal variables is 
P i = E"=i vj W ij with i = 1,... ,n and q ij = v j W i ' n+j - » ! P >+J 
with i,j = l,...,n. We can write 1Z Z in terms of p l , q lJ as 

V' i = 2z(V i - P % H'-=iz 2 (^ V -^ q -- q ^), 

s 

omitting the rather cumbersome linear system for p 1 , q %3 . Similarly as 
in the preceding case, detailed investigation reveals that the operator 
is not invertible. 

Hierarchy of symmetries generated by the operator C is ex- 
pected to contain infinitely many commuting symmetries, as is usual 
with other integrable equations. However, since proof of this fact is 
still pending, we present just first members of the hierarchy. 

As explained in Remark^above, application of the Guthrie opera- 
tor C to a seed symmetry S involves solution of the equation (jlUj) with 
coefficients given by ® with V n \H' %3 replaced by the components of 
the seed symmetry. 



Applying C to the zero symmetry, we infer from Q, ()10|) that Q l 
are constants, henceforth denoted by — c\ Then from we obtain 
an n-dimensional algebra of first-order symmetries 

V i = ^2( V S h Si C S - l^V), 
s 

H ij = tfPj + h'/r ! + h is h sj c s , 

s 

parametrized by c , . . . , c n . By virtue of (J2J) these can be identified 
with the obvious x-translations V 1 = v l s c s , H %3 = ^2 S h l s s c s . These 
were proved by Tenenblat and Winternitz [2Zj to exhaust the algebra 
of Lie symmetries when K ^ 0. 

Application of C to the ^-translations requires solution of system 
fP under substitutions V* = E^v'ft'V - v s h is c!>), H H i = hfd + 
h l - ci + ^2 s h ls h s i c s . With Q % k determined from Q, the potentials Q l 
are just local functions 

The resulting symmetry then has the following components 

s 

+ (2 ^0 ( 3i ^) 2 + Z)(^ S ) 2 + 3^(^) 2 ) ci 
+ 2 J] fr>V + 2 ^ ^ /i st - hf h si )c s 

s St 

+ Y J V s h si (3K(v s ) 2 + J2( hSt ) 2 + 3^(/i' s ) 2 )c s 

s t ( 

and 

s s 

+ hy (sK^f + 3 J^(/i is ) 2 + ^(/i si ) 2 p 

s 

+ tit (3K(vi) 2 + 3 Y,( hjS f + E( /lSJ ) 2 ) c? 
+ 2j2(h%h s i -hi s h*i + h s s ih is y 

s 

+ hiS h Sj ( K (v s ) 2 + J2( hSt ) 2 + 5Z(^) 2 ) cS 

s t t 



9 



(integration constants of (|lUj) are now omitted; they would just gen- 
erate additional x-translations). 

Consider the basis of symmetries V£, IC2 obtained by substitu- 
tion Cj = 5ik- Interestingly enough, expressions V fc *, contain only 
derivatives with respect to x k for each k. In particular, {v l T = Vl, hr = 

} is a system of evolution equations in r and x k , henceforth referred 
to as a generalized pmKdV equation. The motivation is that for n = 2, 
eq. (J2J) becomes the sine-Gordon equation w y i y 2 = — Ksinw under 
identification v 1 = cos ^w, v 2 = sin^w, h 12 = ^w x i, h 21 = — hw x 2, 
y 1 = x 1 + x 2 , y 2 = x 1 — x 2 ;26j- The same identification makes V[ 
into the pmKdV equation wt = w xxx + ^w x (with additional remov- 
able term 3Kw x ). The pmKdV equation is known to coincide with 
the first flow of the hierarchy associated with the sine-Gordon equa- 
tion ^Hl- We conjecture the generalized pmKdV equations and their 
higher analogues to be integrable. 

Hopes are that the hierarchy of generalized pmKdV equations in 
two independent variables possesses a bi-Hamiltonian formulation and 
a classical recursion operator, and the locality of this hierarchy can be 
proved using the results of Sergyeyev j2E] . These and other interesting 
questions will be addressed elsewhere. 

Acknowledgements. Both authors gratefully acknowledge sup- 
port from GACR under project No. 201/04/0538 and from MSMT 
under grant MSM 4781305904. 

References 

[1] Yu. Aminov, On immersions of regions of the n-dimensional 
Lobachevsky space into (2n — l)-dimensional Euclidean space, 
Dokl. Akad. Nauk SSSR 236 (1977) 521-524 (in Russian); En- 
glish transl. it Sov. Math. Dokl. 18 (1977) 1210-1213. 

[2] H. Baran and M. Marvan, A conjecture concerning nonlocal terms 
of recursion operators, this issue. 

[3] J.L. Barbosa, W. Ferreira and K. Tenenblat, Submanifolds of con- 
stant sectional curvature in pseudo-Riemannian manifolds, Ann. 
Global. Anal. Geom. 14 (1996) 381-401. 

[4] R. Beals and K. Tenenblat, An intrinsic generalization for the 
wave and sine-Gordon equations, in: B. Lawson et al., eds., Dif- 
ferential Geometry Pitman Monographs 52 (Longman, 1991) 25- 
46. 



10 



[5] M. Boiti, J. J. -P. Leon and F. Pempinelli, Canonical and non- 
canonical recursion operators in multidimensions, Stud. Appl. 
Math. 78 (1988) 1-19. 

[6] P.T. Campos and K. Tenenblat, Backlund transformations for a 
class of systems of differential equations, Geom. Fund. Anal. 4 
(1994) (3) 270-287. 

[7] J. Cieslinski and Yu.A. Aminov, A geometric interpretation of 
the spectral problem for the generalized sine-Gordon system, J. 
Phys. A: Math. Gen. 34 (2001) L153-L159. 

[8] J. Douglas, Solution of the inverse problem of the calculus of 
variations, Trans. Amer. Math. Soc. 50 (1941) 71-128. 

[9] W. Ferreira, On metrics of constant sectional curvature, Matem. 
Contempordnea 9 (1995) 91-110. 

[10] D. Ferus, and F. Pedit, Isometric immersions of space forms and 
soliton theory. Math. Ann. 305 (1996) 329-342. 

[11] A.S. Fokas and P.M. Santini, Recursion operators and bi- 
Hamiltonian structure in multidimensions. I, Commun. Math. 
Phys. 115 (1988) 375-419. 

[12] A.S. Fokas and P.M. Santini, Recursion operators and bi- 
Hamiltonian structure in multidimensions. II, Commun. Math. 
Phys. 116 (1988) 449-474. 

[13] C.H. Gu and H.S. Hu, Explicit solutions to the intrinsic general- 
ization for the wave and sine-Gordon equations, Lett. Math. Phys. 
29 (1993) 1-11. 

[14] G.A. Guthrie, Recursion operators and non-local symmetries, 
Proc. R. Soc. London A 446 (1994) 107-114. 

[15] N.G. Khor'kova, Zakony sokhraneniya i nelokal'nye simmetrii, 
Matem. Zametki 44 (1988) (1) 134-144. 

[16] I.S. Krasil'shchik and A.M. Vinogradov, eds., Symmetries and 
Conservation Laws for Differential Equations of Mathematical 
Physics, Translations of Mathematical Monographs. Vol.182, 
American Mathematical Society, Providence RI, 1999. 

[17] I.S. Krasilshchik and A.M. Vinogradov, Nonlocal symmetries and 
the theory of coverings: An addendum to A.M. Vinogradov's 
'Local symmetries and conservation laws', Acta Appl. Math. 2 
(1984) 79-96. 

[18] I.S. Krasil'shchik and A.M. Vinogradov, Nonlocal trends in the 
geometry of differential equations: symmetries, conservation 
laws, and Backlund transformations, Acta Appl. Math. 15 (1989) 
161-209. 



11 



[19] S. Kumei, Invariance transformations, invariance group transfor- 
mations, and invariance groups of the sine-Gordon equations, J. 
Math. Phys. 16 (1975) 2461-2468. 

[20] M. Marvan, Another look on recursion operators, in: Dif- 
ferential Geometry and Applications, Proc. Conf. Brno, 
199 5 (Masaryk Univers ity, Brno, 1987) 393-402; ELibEMS 
|http:// www.emis.de/proceedings 

[21] M. Marvan, Some local properties of Backlund relations, Acta 
Appl. Math. 54 (1998) 125. 

[22] M. Marvan, Scalar second-order evolution equations possessing 
an irreducible s^-valued zero-curvature representation, J. Phys. 
A: Math. Gen. 35 (2002) 9431-9439. 

[23] M. Marvan and A. Sergyeyev, Recursion operator for the 
Nizhnik-Veselov-Novikov equation, J. Phys. A: Math. Gen. 36 
(2003) L87-L92. 

[24] P.J. Olver, Evolution equations possessing infinitely many sym- 
metries, J. Math. Phys. 18 (1977) 1212-1215. 

[25] A. Sergyeyev, Why nonlocal recursion operators produce local 
symmetries: new results and applications, J. Phys. A: Math. Gen. 
38 (2005) 3397-3407. 

[26] K. Tenenblat, A note on solutions for the intrinsic generalized 
wave and sine-Gordon equation, J. Math. Anal. Appl. 166 (1992) 
288-301. 

[27] K. Tenenblat and P. Winternitz, On the symmetry groups of the 
intrinsic generalized wave and sine-Gordon equations, J. Math. 
Phys 34 (1993) 3527-3542. 

[28] K. Tenenblat and Chuu Lian Terng, A higher dimension gener- 
alization of the sine-Gordon equation and its Backlund transfor- 
mation, Bull. Amer. Math. Soc. (N.S.) 1 (1979) 589-593. 

[29] C.L. Terng, A higher dimension generalization of the sine-Gordon 
equation and its soliton theory. Ann. of Math. Ill (1980) 491- 
510. 

[30] C.L. Terng, Soliton equations and differential geometry, J. Diff. 
Geom. 45 (1997) 407-445. 

[31] V.E. Zakharov and B.G. Konopelchenko, On the theory of recur- 
sion operator, Commun. Math. Phys. 94 (1984) 483-509. 

[32] Z.X. Zhou, Darboux transformations for the twisted so(p, q) sys- 
tem and local isometric immersions of space forms, Inverse Prob- 
lems 14 (1998) 1353-1370. 



12 



